Magnetic interface states in graphene-based quantum wires 
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The electronic states of a finite-width graphene sheet in the presence of an electrostatic confining 
potential and a perpendicular magnetic field are investigated. The confining potential shifts the 
Landau levels inside the well and creates current-carrying states at or close to the interface with 
the barriers in addition to the edge states caused by the finite width of the sheet. Detailed energy 
spectra are given as a function of the quantum wire parameters. The dependence of the density of 
states on the confinement potential is evaluated for finite and zero magnetic field. 
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I. INTRODUCTION 

The recent production of single layers of stable car- 
bon crystalsiSii has attracted a large interest in their 
fundamental properties and their potential technologi- 
cal applications. The unusual properties of carriers in 
graphene are a consequence of the gapless and approx- 
imately linear electron dispersion at the vicinity of the 
Fermi level at two inequivalent points of the Brillouin 
zone. In the low-energy limit the quasiparticles in these 
systems are described in terms of massless chiral relativis- 
tic fermions governed by the Dirac equation. In particu- 
lar, graphene has been shown to display an unusual quan- 
tum Hall effect ii^Sii, in which the quantum Hall plateaus 
are found in half-integer multiples of 4. This results from 
the fourfold degeneracy of the Landau levels (LL) along 
with the existence of a non-zero Berry phase. The edge 
states in graphene in a perpendicular magnetic field were 
also found to display unusual properties, such as coun- 
terpropagating spin-polarized modest, and are expected 
to play a particularly important role in thin graphene 
structures. 

Another important result concerning single graphene 
layers is the possibility of controlling the electron den- 
sity and Fermi level by a gate voltage. This allows 
the creation of graphene-based quantum structures such 
as potential barriers and quantum wires (QW). Theo- 
retical studies have shown that the relativistic behav- 
ior of quasiparticles in graphene allow the observation 
of effects such as the Klein paradox, which is the per- 
fect transmission of relativistic particles across poten- 
tial barriers, as well as a direction-dependent tunneling 
through barriers2*i2*ii*i2iiii. In addition, recent exper- 
imental work has demonstrated electronic confinement 
in patterned graphene structures created by standard 
lithography methods^. 

In this paper we study the interplay of an electrostatic 
potential barrier and an external perpendicular magnetic 
field on a graphene QW and find that propagating states 
exist at the interface with the potential barriers in addi- 
tion to the ordinary edge states. Further, for sufficiently 
wide potential wells two distinct sets of LL arise due to 



the energy shift caused by the presence of barriers. 

This paper is organized as follows. In Sec. II we 
present the model and formalism and in Sec. Ill nu- 
merical results. We conclude with remarks in Sec. IV. 



II. MODEL AND FORMALISM 

The crystal structure of undoped, defect-free graphene 
layers is that of a honeycomb lattice of covalent-bond car- 
bon atoms. To each carbon atom corresponds a valence 
electron and the structure can be described as composed 
of two sublattices, labelled A and B. The low-energy ex- 
citations of the system at the vicinity of the K point and 
in the presence of both an electrostatic potential U and 
a uniform magnetic field B perpendicular to the plane 
of the graphene sheet are described, in the continuum 
approximation, by the 2D Dirac equation 

{v F [a-p- eA]+mvpa z }^ = (E - U)%, (1) 

where the pseudospin matrix a has components given by 
Pauli's matrices; p = (p x ,p y ) is the momentum opera- 
tor. The "speed of light" of the system is vp, the Fermi 
velocity (vp w 1 X 10 6 m/s), and A is the vector po- 
tential. The eigenstates of Eq. (1) are represented by 
two-component spinors \& = [tpA > V , -b] T , where ipA and 
vjb are the envelope functions associated with the proba- 
bility amplitudes at the respective sublattice sites of the 
graphene sheet. The term cx mv F introduces an energy 
gap, which may represent e.g. the effect of spin-orbit 
coupling. 

We now consider a narrow graphene layer, of width 
W, in the presence of a one-dimensional (ID) potential 
U = U{x) and a perpendicular magnetic field B. This 
allows us to write the solutions for the spinor compo- 
nents in the form ipA(x,y) = (f>A(x)e %kyV and ipB(x,y) = 
i4>B{x)e lkyV because of translational invariance along the 
y direction and the particular choice of Landau's gauge 
A = (0, Bx,0). The resulting equations for 4>a{x) and 
4>b{x) are 
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+ {k y -eBx)4>B = [E-U(x)-mv 2 f ]<l> A . (2) 



These equations can be decoupled and, by setting £ = 
/3 1 / 2 (x — ky/P), where /3 = £ B 2 — eB/h is the inverse 
magnetic length squared, the result is 
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[d(t> A /dt + £<j> A ] = 0, (3) 



[*M-W=0, (4) 



where fi = [(e - u) 2 - A 2 ]//3, w = (tiv F )- 1 U, e = 
(Hvf)~ 1 E, A = mvp/h and the prime denotes deriva- 
tive with respect to £. For a constant potential U = Uo, 
Eqs. (2) and (3) have well-known solutions in terms of 
Hermitc polynomials and the spectrum is 



E = ±hv FV / 2nP + A 2 + Uo, 



(5) 



where n is an integer. This contrasts significantly with 
the nonrelativistic spectrum E = huj c (n + 1/2) + Uo- 



III. RESULTS 

A. Dispersion relation and density of states 

First we consider the effect of a steep potential well 
U(x), with a characteristic width L in a graphene strip 
of width W . In this case the derivatives of the potential 
are strongly localized functions that have non-zero values 
only at the vicinity of the barrier interfaces. The solu- 
tions then depend on the width W and the strength of 
the magnetic field B through the ratio £b/L. Let us con- 
sider initially the case 2£b/L « 0.5, which corresponds 
to classical orbits that fit inside the QW. We can assume 
solutions of the form 4>c = /c(£) e_ ^ 2 ' C — A,B. For a 
constant potential fc are the well-known Hermite poly- 
nomials. Thus, for sufficiently strong B and small k y the 
spinor functions quickly decay with £ and the solutions 
are, to a good approximation, localized inside the well. 
Therefore, the energy spectrum is dispersionless and for 
small n is that of the LL for U = 0. For larger values 
of k y the center of the solutions £ 2 B k y is shifted towards 
the barrier regions. For the lowest LL (n = 0) one can 
estimate the limits of the central dispersionless region by 
setting exp(— £ 2 /2) w 0.1, or £ ps 2, since in this case the 
amplitude of the wave function inside the barriers is neg- 
ligible. This gives -1 + U B /L < 2k y £%/L < 1 - U B /L. 
For non-zero values of n, the central dispersionless re- 
gion of the spectrum is expected to be narrower, since 
the spinor functions are less localized and can have a 
larger magnitude within the barrier regions. 



For sufficiently large values of the momentum along 
the y direction, a similar argument shows that the en- 
ergy levels may be accurately approximated by the LL 
shifted by U in the regions 2k y £ 2 B /L < -1 - 4£ B /L 
and 2kyl 2 B jL > 1 + 4£b/L. For intermediate values 
of k y one expects dispersive solutions. These solutions 
can be described as interface states, in the sense that 
they are localized electronic states in the x direction, 
that propagate along the interfaces with the potential 
barriers, in analogy with the edge states of a 2D elec- 
tron gas in a magnetic field, but with the fundamental 
difference that in the present case the spinor functions 
are non-negligible both inside the QW and in the barri- 
ers. These are current-carrying states and, for sufficiently 
smooth potentials, they should not depend on the micro- 
scopic structure of the graphene sheet. For even larger 
values of k y the spinor functions may be shifted toward 
the edges of the sample and give rise to edge states, which 
have been shown to depend on the shape of the graphene 
edgesi^. 

We have considered two specific types of potential 
wells, as shown in Fig. 1. The solid line illustrates the 
step potential given by 

U{x) = ^{t a nh[(-x-^)/S}+t & nh[(x-^)/5}+2}, (6) 

where 5 denotes the thickness of the interface. The 
dashed line in Fig. 1 refers to a parabolic potential 

( U (2x/L) 2 \x\ < L/2, 
U{x) = { (7) 
[ Uo \x\ > L/2. 

Figure 2 shows numerical results for the energy spectrum 
of a QW with the potential given by Eq. (6). The re- 
sults are for L = 200 nm and an interface width of 10 
nm, which is much larger than the lattice parameter of 
graphene a sa 0.14 nm, B = 2 T and Uq = 50 meV. In 
this case we have 2£b/L « 0.18. The vertical dashed lines 
delimit the range 1 - Al B /L < \2k y £ 2 B /L\ < 1 + U B /L. 
In contrast with the conventional edge states, the bar- 
rier interface does not cause a splitting of the disper- 
sion branches. This is a consequence of the difference in 
boundary conditions at the potential interfaces and the 
edges of the graphene sheet, i.e., at the potential step 
the wave function is finite and continuous and the associ- 
ated probability density can be significant for \x\ > L/2, 
whereas at the edges of the sample, as in the case of a 
graphene sheet with an armchair termination, the wave 
function was assumed to vanish. 

For comparison, Fig. 3 shows the spectrum of a 
parabolic confining potential given by Eq. (7). This po- 
tential is non-zero at every point, except at x — 0, and is 
the reason why the states are dispersive for all values of 
k y satisfying \2k y £ B /L\ < 1. This result becomes similar 
to that of the tangent hyperbolic potential for large S. 

For weaker fields, such that 2£b/L > 0.5, the former 
picture breaks down and the magnetic field acts as a per- 
turbation to the zero-field case 10 . For B — the electron 
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FIG. 3: The same as Fig. 2 but for a parabolic potential given 
by Eq. (7), with U = 100 meV, L = 200 nm, B = 2 T. 
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FIG. 2: Energy spectrum of a graphene QW for B = 2T and 
a potential given by Eq. (6), with Uo = 50 meV, L — 200 nm 



states inside the QW that propagate perpendicularly to 
the barrier interfaces are transmitted without reflection 
(Klein tunnclling)2ii2*ii*i£. This counterintuitive behav- 
ior results from the absence of a gap in the spectrum and 
from the chiral nature of the quasiparticles in graphene. 
However, recently it has been demonstrated that non- 
zero values of momentum along y allow the existence of 
confined electron states in a QW— . For large values of 



E = hv F [{£w/L) 



k 



211/2 



(8) 



where £ is an integer. As a finite magnetic field is intro- 
duced, one can expect a modification of these states. In 
particular, one expects the existence of localized states 
at k y = 0. This situation is observed in Fig. 4, where the 
energy spectrum is plotted as a function of wave vector, 
for the potential given by Eq. (6), with Uo = 100 meV, 
L = 100 nm, W = 800 nm and B = 0.6 T. This case 
corresponds to 2£b/L ~ 0.66. The figure shows the exis- 
tence of dispersive states for a wide range of wavevectors, 
as well as dispersionless branches that correspond to the 
shifted LL at the barriers. 

For small energies and wavevectors, the results for B = 
show the existence of hole states in the barriers that 
either propagate (Klein paradox) or tunnel through the 
well. For a finite field, in a semiclassical description, 
the trajectories of the holes would be deflected by the 
magnetic field and would be confined into closed orbits 
that cross the QW. A similar behavior is thus obtained 
for energies close to zero, as shown in Fig. 5. The figure 
shows the wave functions (left, panels (a) and (c)) and 
the respective probability densities (right, panels (b) and 
(d)) for two states with k y = in a QW with the same 
parameters as in Fig. 4 and two energies: E = 2.03 meV 
(upper panels) and E = (lower panels). The higher- 
energy state describe holes that cross the well region via 
electron-hole conversion and thus show a maximum in the 
probability density inside the well, whereas the result for 
zero energy corresponds to holes that tunnel across the 
well region and are confined by the magnetic field. For 
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FIG. 4: The same as Fig. 2 but now with U = 100 meV, 
L = 100 nm and a smaller magnetic field B = 0.6 T. 
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FIG. 5: Color online. Wave functions (left panels, <j>A'- solid, 
4>b'- dashed-dotted) and the respective probability densities 
(right panels) for two low-energy states of Fig. 4. Upper 
panels: E = 2.03 meV, lower panels: E — meV. 



larger energy eigenstates and k y = the wave functions 
show an oscillatory behavior inside the potential well and 
quickly decay in the barriers. Therefore, these states 
possess a dominant electron-like character. 

For finite values of k y the spectrum displays dispersive 
branches with either positive or negative slopes. The 
wave functions in these cases have non-negligible ampli- 
tudes, respectively, inside the potential well and in the 
barriers, and thus are associated with propagating elec- 
trons and holes. On the other hand, the states associated 
with flat energy branches have finite amplitudes inside 
the barriers. The propagating states are found to inter- 
act with the dispersionless states, with the appearance of 
anticrossings. 

The right part of Fig. 6 shows comparative plots of 4>a 
for electron sates at the proximity of a specific anticross- 
ing for the situation corresponding to Fig. 4 and blown 
up on the left part of the figure. The results for points 
(a) and (b) (k y 2£ B /L = 2.19) indicate that close to the 
anticrossings, the electron states are a superposition of 
oscillatory states (inside the QW) and non-propagating 
states (at the barriers), whereas away from the anticross- 
ings (k y 2£ B /L = 2.59) the wavefunctions match the re- 
sults for either confined states (c) in the QW and LL (d) 
in the barriers. The vertical lines (dotted) delimit the 
QW region. 

In the calculation of several physical quantities the 
density of states (DOS) is needed, which is defined as 
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FIG. 6: Left: A zoom of the electron dispersion branches of 
Fig. 4 at the vicinity of an anticrossing; the parameters are 
given in the text. Right (a-d): Spinor component <f>A for some 
particular states specifed in the text. 



(9) 



fc x ,fc„ n 



where S = WL is the total area of the sample and n 
labels the different energy branches. The summation over 
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k x is replaced by an integral, and we obtain 

where we introduced broadening of the energy levels by 
replacing the 8 functions with Lorentzians of constant 
width T: Eq = \/2ehB denotes the characteristic energy 
scale of the system. 

Figure 7 shows numerical results for the DOS for the 
hyperbolic tangent QW, with Uq — 50 meV, L = 200 nm, 
and for different values of the external magnetic field: 
B = 0.5 (dashed), 1.0 (dotted), 1.5 (dash-dotted) and 
2.0 T (solid). In all cases we used T/E = 28. The 
figure shows that the presence of the barriers shifts the 
DOS peaks by Uq from the zero potential results. In 
contrast with a conventional 2D electron gas, all DOS 
results for graphene show a pronounced peak at E = 50 
meV whereas the remaining peaks arc shifted according 
to the different values of B. This is a consequence of the 
fact that in graphene, the energy of the LL with index 
n = is independent of B. In addition, for each value of 
B extra peaks result from the LL inside the QW. Due to 
the square root dependence of the spectrum, see Eq. (5), 
the distance between the peaks decreases as the energy 
increases. Therefore, the influence of the LL in the QW 
becomes more evident for energies closer to zero. The 
inset contrasts the DOS of a non-relativistic electron gas 
(green) with the DOS of electrons in graphene (solid), for 
Uq = 0. 

The dependence of the DOS on the barrier height Uq 
is shown in Fig. 8, where L — 200 nm and B = 2 T, for 
different values of Uq, namely Uq — (solid line), Uq = 25 
meV (dashed), Uq = 50 meV (dotted) and Uq = 100 
meV (dot-dashed). As the potential increases, there is a 
clear shift of the peaks from the values given by Eq. (5) 
with Uq = 0, with the n = peak now shifted to the 
value of Uq. Also evident is the presence of states inside 
the QW, indicated by the existence of additional peaks, 
with positions that are independent of the potential step. 
These results indicate that, due to the specific nature 
of the LL spectrum in graphene, one can increase the 
DOS at a particular energy by a suitable change of Uq so 
that one LL outside the QW is matched with another LL 
inside the well. This condition (for A = 0) is expressed 
as 



U = hv F (^7p - ^0), (11) 

where n and n' are integers. This condition is approached 
in Fig. 8 in the result for Uq = 25 meV, for the peak at 
E = 84 meV, with n! = 4 and n = 2. Figure 9 shows 
the DOS as a function of the external magnetic field for 
E = 8 meV, for Uq = (solid line) and Uq = 50 meV 
(dashed line). The shift of the LL brought about by the 
potential barriers causes the appearance of several DOS 
peaks that are absent from the uniform system. 
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FIG. 7: Color online. Density of states for a QW on graphene 
with U = 50 meV, L = 200 nm. B = 0.5 T (dashed), 
B = 1.0 T (dotted), B = 1.5 T (dash-dotted) and B = 2.0 
T (solid). Inset: comparison between the DOS of a non- 
relativistic electron gas (green) and the DOS of electrons in 
graphene (solid), for Uo = 



B. Energy spectrum versus QW parameters 

Figure 10 shows the spectrum, as a function of the 
magnetic field, for two values of the wave vector, k y — 
(a) and k y — 0.1 nm -1 (b), for the QW potential of Eq. 
(6) with Uq = 100 meV and L = 100 nm. In both cases, 
the results show that the presence of the potential bar- 
riers introduces a significant modification of the energy 
eigenstates, in comparison with the results for Uq = 
(in which case the eigenvalues are proportional to the 
square root of the external magnetic field). Figure 10(b) 
shows the presence of the quantized confined states as B 
tends to 0. These discrete states are initially very weakly 
dependent on the external field. As the field increases, 
these states interact with the non-propagating states, as 
evidenced by the presence of anticrossings. At small fields 
the figure shows split branches, which are edge states that 
arise due to the fact that the spinor amplitudes are set to 
zero at the edges of the sample, whose width is W = 800 
nm. This corresponds to the boundary conditions of a 
graphene sheet with an armchair edge. 

The left panel in Fig. 11 shows the dependence of 
the energy eigenvalues on the QW width, for k y = 0, 
B = 2 T and Uq = 100 meV. The results show that the 
higher-energy levels are strongly modified even for wide 
wells, whereas the lower-lying states can be accurately 
described by Eq. (5) with Uq = for a relatively thin 
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FIG. 8: Color online. Density of states for a graphene QW 
with B = 2 T and L = 200 nm. The solid, dashed, dotted, 
and dash-dotted curves correspond, respectively, to Uo = 0, 
U = 25 meV , Uo = 50 meV, and Uo = 100 meV. The last 
three curves are shifted up by 1.25 for clarity. 
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FIG. 10: Energy spectrum of a graphene QW as a function of 
the magnetic field B, with Uo = 100 meV, L = 100 nm, (a) 



k y — 0, and (b) k y =0.1 nm . 
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FIG. 9: Color online. Density of states as a function of the 
magnetic field for E = 8 meV, in the absence of a confining 
potential (black solid curve) and a graphene QW with Uo = 50 
meV (blue dashed curve), L = 200 nm in both cases. 
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FIG. 11: Energy spectrum of a graphene QW for B — 2 
T and k y = 0. In the left panel the spectrum is plotted vs 
the width L, for Uo = 100 meV, and the right panel vs the 
potential height Uo for L = 100 nm. 
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QW. The figure allows us to distinguish two regions in 
the spectrum, which arise due to the non-uniform distri- 
bution of the LL in graphcnc. One region corresponds to 
the lower-energy states, which arise due to the strong in- 
teraction between the low-energy states in the well with 
the shifted negative-energy states in the barriers. The 
degeneracy of these states is lifted for small L and the 
spectrum shows an approximately linear dependence on 
the QW width. In particular, the interaction with the 
lower-energy states modifies the E = LL for L ss 80 
nm, which is equivalent to 2£ B /L ss 0.45. The other re- 
gion corresponds to higher-energy states, which remain 
degenerate but are strongly shifted for small L. These 
states are weakly dependent on the well width for larger 
values of L, e.g., for L > 120 nm. 

The dependence of the energy spectrum on the po- 
tential barrier height Uo is shown in the right panel of 
Fig. 11 for L = 100 nm, B = 2 T, and k y = 0. As 
Uo increases, the figure again shows the appearance of 
two sets of states: one comprises states that are weakly 
dependent on Uo and the other states that show a signifi- 
cant dependence on Uo- This is caused by the hole states 
in the barriers, whose energies are shifted by the poten- 
tial. As Uo increases, the LL in the well region interact 
with the set of shifted of LL in the barriers, causing the 
appearance of additional electron states at energies close 
to zero. 



a graphene QW in the presence of a perpendicular mag- 
netic field. We found a shift of the Landau levels caused 
by this potential which may be observable by its effect 
on the quantum Hall steps in the presence of gate volt- 
ages. For steep potential barriers, there is a clear distinc- 
tion between a low magnetic-field regime, characterized 
by the existence of dispersive confined states with small 
wave vectors inside the QW as well as non-propagating 
states in the barriers, and a higher magnetic-field regime, 
in which there are non-propagating states inside the well 
and outside the QW, together with propagating states at 
the interfaces of the potential barriers. These interface 
states may cross the Fermi level and, together with the 
conventional edge states that arise due to the finite size 
of the sample, can contribute, e. g., to the conductivity 
of the system. The modification of the LL spectrum in 
the QW is also evident in the shift of the peaks in the 
DOS in comparison with the results for zero confining 
potential, along with the appearance of additional peaks 
caused by the LL inside the QW. 
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